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We use the density functional method to examine the properties of the nonuni- 
form (two-phase) fluid of two-level atoms, a part of which is excited. From the 
analysis of the equation of state of a gas of two-level atoms, a part of which is 
excited, the following density functional of the grand thermodynamical potential 
emerges 
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(here £lcs[p(r)] is the Carnahan-Starling term, a is the atom radius, v = |7rcr 3 , 
ci is the concentration of excited atoms, cq + c\ — 1, E\—Eq is the excitation 
energy and a is the dimensionless parameter which characterizes the atom). We 
use this expression to calculate the nucleation barrier for vapor-to-liquid phase 
transition in the presence of excited atoms. 
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The studies of the equilibrium properties of a gas of identical atoms, a part 
of which is in an excited electronic state, attract attention during last forty years [ 
d El El HI El El 13 El HE CDj • Such atoms may appear owing to the electromagnetic irra- 
diation with the frequency, which corresponds to the excitation energy of the atom. 
Since the life-time of the excited state is essentially larger than the time required 
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for establishing of equilibrium over translational degrees of freedom, the system 
should exhibit equilibrium properties at a given (nonequilibrium) concentration of 
excited atoms. Moreover, owing to new effective long-range interatomic interactions 
- the resonance dipole-dipole interactions - one may expect essential changes of var- 
ious equilibrium characteristics due to the presence of excited atoms. The analysis 
performed within the frames of the cluster expansion method confirmed these expec- 
tations (see Refs. [0IE1 and references therein). On the other hand, there are only 
a few papers which may be related to experimental observations of the theoretically 
investigated features of the gas with excited atoms. We should mention here the 
papers on the influence of irradiation on the condensation of iodine and anthracene 
vapor [ HU H21 and on the photonucleation [ H31 [T^l [THl [THl [TTl HH1 HH1 1201 12H , in 
particular, in vapors of mercury and cesium. The latter studies reported the quan- 
titative results of resonance irradiation influence on the nucleation rate. The main 
conclusion of these studies is as follows: the resonance irradiation of the nucleation 
zone leads to a sharp increase of the nucleation rate. The obtained experimentally 
results, to our best knowledge, have not been explained until now. One may try to 
interpret these data basing on the theory of nucleation in a supersaturated vapor 
which contains excited atoms. Such atoms may appear as a result of the resonance 
irradiation. 

In what follows we present preliminary results about the nucleation phenomena 
in a fluid of two-level atoms, a part of which is excited (see also [122])- In our study 
we use the density functional approach developed by D. W. Oxtoby with coworkers [ 
23J. This method permits to obtain the nucleation rate basing on the first principles. 
The nucleation rate J is connected with the nucleation barrier A, J = J exp ( _ ^) • 
The nucleation barrier can be calculated within the frames of the classical nucleation 
theory which relies on the capillarity approximation [ 123] 

A d _ 167T / 7 y 1 

kT 3 KkTj pf\n 2 s' { ' 

Here 7 is the surface tension of the vapor-liquid interface, pi is the density of liquid, 
s = 2- is the supersaturation, p is the actual pressure of the supersaturated vapor 
and po is the equilibrium pressure. To compute the nucleation rate according to Eq. 
(PQ) one has at first to construct the vapor-liquid phase diagram determining the 
equilibrium pressure po and the liquid density pi at the temperature T and then 
to compute the surface tension 7 at this temperature. This calculation can be done 
within the density functional approach considering the planar vapor-liquid interface. 
Alternative approach for calculation of the nucleation barrier which does not use the 
key assumption in the classical nucleation theory - the capillarity approximation - 
was suggested by D. W. Oxtoby According to this scheme one has to consider 
a metastable vapor in a spherical vessel of the radius 1Z, assume appearance of a 
spherical liquid droplet in the center of the vessel and analyze the density profile and 
the value of the grand thermodynamical potential of such a two-phase fluid with the 
spherical vapor-liquid interface. The value of the grand thermodynamical potential 
of such a metastable fluid Q(T,p, V) permits to calculate the nucleation barrier via 
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the equation 

A = Q(T,fi,V)- (-ptnT&y (2) 

In our study we use both schemes for calculation of the vapor-to-liquid nucleation 
barrier in a system of two-level atoms, a part of which is excited. 

To perform the theoretical analysis of the vapor-to-liquid nucleation in the pres- 
ence of the excited atoms we need an appropriate density functional of the grand 
thermodynamical potential. We assume the grand t her mo dynamical potential to be 
a functional of the density with the form 



Q[p(r)] = kT J dr lP (n) (ln(A 3 p(n 
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where 



a(d, T) = ^v(Ei - E ) (l - 2cj + 2(1 - c^c^^) . (4) 

Here A is the thermal de Broglie wavelength of the atom, v = f^c 3 , 0" is the radius of 
the atom, c\ is the concentration of the excited atoms, a = %J° E is the dimensionless 
parameter which characterizes the two-level atom (in what follows we set a = 1), 
Ei — Eq is the excitation energy, d is the value of the transitional electrical dipole 
moment between the ground and excited states. The equilibrium density minimizes 
the grand thermodynamical potential fl[/?(r)], i.e. is the solution of the following 
integral equation 

kTln (AVn)) + tf Mri1 " 9l ' 2p2(ri) + 3 3 " V(r0 

nvMcuT) f dr2 ^M .^o. (5) 
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Substituting the equilibrium density into Eq. Q one obtains the value of the grand 
thermodynamical potential of the system Q(T, fi, V). The adopted density functional 
of the grand thermodynamical potential (J3J) is consistent with the virial state equa- 
tion obtained earlier [ |T] . It takes into account the short-range interaction within 
the Carnahan-Starling local approximation (the first term in the r.h.s. of Eq. ©) 
neglecting the difference of the atom radii in the ground and excited states. More- 
over, it takes into account the long-range interactions, in particular, the resonance 
dipole-dipole interactions, within the mean-field approximation (the second term in 
the r.h.s. of Eq. (jSJ)). Note, that the coefficient a(ci,T) (@J depends on temperature 
only when c\ deviates from zero. More sophisticated density functionals are available 
but they have not been employed in the present study. 
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It is convenient to introduce the dimensionless units of energy, temperature, 
chemical potential, length, volume, density, pressure, surface tension etc renormal- 
izing these quantities as follows: E -> T -> p -> ^4^, r ~+ 

y -> p -> up, p -> 7 -> etc ' res P ectivel y- As a result E q s - ©, O, 

become 

3T r . / /A 3 A -l + 6p(r 1 )-4p 2 (r0 
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«(ci,T) = — I 1 - 2ci + I, (7) 



rin f^ (ri) V T Mr 1 )-9p-(r 1 ) + 3p3(r 1 ; 



u / (I-P(ri)r 
lMCl ' T) ^ dr 2 ^^-p = 0. (8) 
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Moreover, we assume for concreteness in Eqs. (@J) and (|7|) a = 1. We also set without 
loss of generality — = 1. 

We start with the phase diagram of the system. For this purpose we assume the 
constancy of the density, p(r) = p, that immediately yields instead of Eqs. JBJ) and 
(JHJ) the following expressions 

= Tp^ (hip + ~ l + 6p ~ 4p2 ^) -a( Cl ,T)pV-ppy (9) 
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and 



Tlnp + T 8p 9p2+ 3 3p3 -2Q( Cl ,r)p-p = 0. (10) 

Solving Eq. (fTUjl with respect to p and substituting this density into Eq. (jHJ) one gets 
the value of the grand thermodynamical potential fi(T, p, V). One can also eliminate 
with the help of Eq. (jlUJl the chemical potential p from Eq. (JSJ) getting as a result 
the equation of state 

n(T,p,V) _p_ l + p + p 2 -p 3 a(ci,T) 2 

Ty t p (i - p) 3 r p ' 1 j 

Eq. (jllj) agrees with the virial state equation of a gas of two-level atoms, a part 
of which is excited, obtained earlier [Q] (see also Refs. [0IH])- Eq- m ay have 
more than one solution which yield the same value of the grand thermodynamical 
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potential f2(T, p, V). Indeed, for a given temperature T let us fix the value of the 
grand thermodynamical potential — y=p and solve Eq. (fTT]) with respect to p. At 
high temperatures (above the critical temperature T c ) one finds only one solution 
p which corresponds to a certain value of the chemical potential p in Eq. (jlOj) . At 
low temperatures (below the critical temperature T c ) one finds several solutions p 
with the corresponding values of chemical potential p which follow from Eq. (jl(J|) . 
Varying the value of the grand thermodynamical potential — U = p one finds such 
two densities p v and pi > p v which yield the same value of p. The quantities T, 
P = Po, Pv, Pi, P = Po correspond to the points on the phase diagram where the 
two phases, liquid and vapor, coexist (see Fig. [TJ). In Fig. [T]we display the isotherms 




Figure 1. The isotherms p vs - at T = 0.6T C (0) (left panel) and T = 0.8T C (0) 
(right panel). Bold curves correspond to c\ = 0, thin curves correspond to 
ci = 0.00006 (= 0.006%). The left (right) endpoint of the horizontal part of 
the isotherm (corresponding to po) gives the liquid density pi (the vapor density 

Pv)- 



which correspond to the temperatures T = 0.6T C (0) ~ 0.00176866 (left panel) and 
T = 0.8T C (0) « 0.00235822 (right panel) (here T c (0) denotes the critical temperature 
T c without excited atoms, i.e. when c\ = 0) for two concentrations of excited atoms, 
Ci = (bold curves) and C\ = 0.00006 (= 0.006%) (thin curves). Considering at first 
the case C\ = at T = 0.8T C (0) we find that the equilibrium values of the pressure, 
the chemical potential, the liquid density, and the vapor density are po ~ 0.00004118, 
~ -0.00996108, pi w 0.30719568, and p v w 0.02172324, respectively. Assume 
further that in a system the concentration of excited atoms becomes c\ = 0.00006. 
For such a fluid the equilibrium values of the pressure is p Q ~ 0.00003179 and 
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the vapor with the pressure ~ 0.00004118 becomes metastable with the value of 
the supersaturation parameter s ~ 1.29559738. Moreover, the equilibrium values 
of the chemical potential, the liquid density, and the vapor density of the fluid 
with ci = 0.00006 at T = 0.8T C (0) are // w -0.01049448, pi « 0.32745048, and 
p v ~ 0.01596817, respectively. 

To calculate the vapor-to-liquid nucleation barrier according to Eq. (JTJ) one has 
to find the surface tension 7. Analyzing the density profile for a planar vapor-liquid 
interface (for this purpose we consider a two-phase system in a cylinder of the radius 
K and the height C) at T = 0.8T C (0) and c x = 0.00006 and estimating Q(T, fa, V) 
we find according to the relation ■ynTl 2 = tt(T, fa, V) — (—p n7l 2 C) the value of the 
surface tension 7 = 0.00051195. As a result one immediately gets the value of the 
vapor-to-liquid nucleation barrier ^ « 68.3294 (see Fig. |2l dash-dotted curve 3). 
Obviously, since ^ becomes now finite (and decreases as C\ increases) the nucleation 




0.00006 0.00008 



Figure 2. The dependence of the vapor-to-liquid nucleation barrier on the con- 
centration of excited atoms c\ at two temperatures T = 0.6T C (0) (curves 1 and 
2) and T = 0.8T C (0) (curves 3 and 4). The curves 1 and 3 were obtained using 
Eq. (0), the curves 2 and 4 were obtained using Eq. (J2J). 

of liquid from vapor becomes now possible. 

On the other hand, we can calculate the nucleation barrier on the basis of Eq. 
(J2J). First we estimate the Thompson radius r* = kT 2 p ^ ln - at T = 0.8T C (0) and 
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Cl = 0.00006. We obtain £ « 6.614283. We note that r* is rather small that 
may be a reason to go beyond the classical nucleation theory since the capillarity 
approximation cannot be justified for such small droplets. Next we calculate the 
chemical potential for the supersaturated vapor with excited atoms according to 
Eq. (JTUJ) with p « 0.02227178 (this value of density follows from Eq. (jnj) for c x = 
0.00006, T = 0.8T C (0) and p w 0.00004118), // = -0.01049448 + 0.00049852, and 
analyze the density profile of a spherical droplet in the supersaturated vapor seeking 
for a "stable" value of the grand thermodynamical potential which plays a role of 
ft(T,//, V) in Eq. © (for details see [E2]). We find | « 67.12 that agrees with the 
value obtained within the frames of the classical nucleation theory. The described 
calculations have to be repeated for other values of concentration c\. Moreover, we 
perform such calculations for several values of temperature. Some of our findings are 
collected in Fig. |21 

The main conclusion which can be read off from Fig. |21 is as follows: the vapor 
after appearance of excited atoms becomes metastable with s > 1 and the nucleation 
barrier for vapor-to-liquid phase transition becomes essentially diminished. This out- 
come agrees with a naive expectation that the long-range resonance dipole-dipole 
interactions should act in favor of liquid formation in vapor. Although the present 
consideration permits to obtain the nucleation rate which can be measured experi- 
mentally much more work is required to compare theory and experiment. First, we 
have to analyze in detail the results for nucleation rates [ H3 EH] obtained using 
the upward thermal diffusion cloud chamber setup [E3|- Second, we should bare in 
mind that the fluids whose photonucleation has been studied have more complicated 
particle structure and interparticle interactions. Comparison with experiment can 
therefore be only qualitative at present, and in this respect our results are consistent 
with the data reported in Refs. [ H3J [TE] . 

One of the authors (O. D.) is grateful to the DAAD for the support of his 
visit to Philipps-Universitat Marburg in the autumn of 1995. He wishes to thank 
Dr. Hermann Uchtmann for kind hospitality and many stimulating conversations. 
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